AD-A053  574 


UNCLASSIFIED 


MISSISSIPPI  UNIV  UNIVERSITY  DEPT  OF  ELECTRICAL  EN6IN— ETC  F/S  20/14 
ANGULAR  VARIATION  OF  CURRENT  ON  MODERATELY  THIN  CYLINDRICAL  SCA— ETC(U) 
SEP  76  B M DUFF  AF0SR-75*2862 

AFOSR-TR-78-0769  NL 


AD  A053574 


^iFOSR-TR-  ? R - 0 7 6 9 


(3^ 


! L I 


MGULAR  VARIATION  OF  gpRRENT  ON 
MODE^TELY  THIN  CYLINDRICAL  SCATTERERS, 


( 1/-^ ■. 

T^'  1, B^M./Duff' 

Department  of  Electrical  Engineering 
University  of  Mississippi 
University,  Mississippi  38677 


(i'V 


D D C 

FfTapnnnc 

MAY  5 1976 


/ , j r ■ ■ ' ^ ^ ''  'W 


Abstract 


The  surface  current  density  induced  on  finite  length  cylindrical  scat- 

terers,  of  moderately  thin  radius,  has  been  investigated  experimentally. 

The  results  of  this  investigation  clearly  show  an  angular  variation  of  the 

axially  directed  current  density  on  cylinders  of  electrical  radius  ka  = 0.04. 

The  data  has  been  analyzed  to  express  the  current  density  as  a Fourier  series 

I ’ r 

in  the  angle  variable  ip.  The  results  of  this  analysis  show  that  the  angularly 
Independent  component  is  essentially  in  agreement  with  solutions  based  on  a 
thin  wire  approximation,  while  the  cos(n()))  terms  are  essentially  constant 
with  respect  to  z.  Comparisons  are  then  made  with  the  known  solutions  for 


an  infinitely  long  cylindrical  scatterer.. 
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(1)  Introduction 

Conducting  structures  composed  of  electrically  thin  cylinders 
(wires)  form  an  important  class  of  problems  in  electromagnetics.  When 
excited  by  a source  within  the  structure,  these  geometries  comprise  the 
linear  antenna  problem.  When  excited  by  sources  exterior  to  the  struc- 
ture, they  are  treated  as  scatterers.  The  thin  wire  scatterer  has  often 
been  used  to  obtain  a first  approximation  to  more  complicated  structures, 
for  example,  the  modeling  of  an  aircraft  with  the  wings,  fuselage,  tail, 
and  so  forth  replaced  by  sections  of  electrically  thin  conducting 
cylinders.  In  connection  with  these  theoretical  solutions,  experimental 
measurements  of  the  current  and  charge  distribution  have  been  made  [1,2]. 
The  importance  of  the  electrically  thin  cylinder  arises  from  simplifica- 
tions which  may  be  made  in  analysis  of  the  structure.  These  simplifications 
are  generally  lumped  under  the  term  "thin  wire  approximations".  It  is 
always  of  some  importance  in  both  theoretical  and  experimental  investi- 
gations to  evaluate  the  effect  of  the  thin  wire  approximation  and  the 
restrictions  on  electrical  radius  which  are  necessary  for  required 
accuracy.  In  the  case  of  driven  antenna  problems,  the  thin  wire  limit 
is  often  taken  as  ka-0.01  where  k is  propagation  constant,  a is  the 
radius  of  the  cylinder.  Experimental  measurements  on  driven  antenna 
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structures  with  ka=0.04  have  been  made  and  show  good  agreement  with 
theoretical  solutions  based  on  thin  wire  approximations  [1].  An  assump- 
tion which  is  common  to  all  thin  wire  approximations  and  the  one  with 
which  the  present  work  is  concerned  is  that  the  current  distribution  in- 
duced on  the  structure  is  independent  of  the  cylindrical  angle  variable 
(p.  As  will  be  shown  this  assumption  is  strictly  valid,  for  cylinders  of 
nonzero  radius,  only  when  the  excitation  is  Itself  Independent  of  (j).  The 
image  plane  monopole  driven  by  a TEM  mode  coaxial  line  Is  an  example  of 
a <|)  independent  excitation.  It  has  been  recognized,  however,  that  in  a 
scattering  problem  for  which  the  excitation  is  an  Incident  electromagnetic 
wave  the  assumption  of  angular  independence  of  the  current  is  not  strictly 
valid. (For  a discussion  of  this  effect  related  to  the  EMP  problem,  see 
reference  3.) 

The  necessity  for  a (()  dependence  of  the  current  distribution  induced 
on  a conducting  cylinder  by  an  Incident  plane  wave  may  be  seen  from  an 
examination  of  the  boundary  conditions  on  the  H field.  Referring  to 
Fig.  1,  the  Important  boundary  condition  is 

n X ii  = j (1.1) 

8 

where 

n = the  unit  normal  to  the  surface 

H = the  total  H field 

J = the  induced  surface  current 
s 
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Expressing  the  total  H field  as  a combination  of  incident  and  scattered 
fields,  we  obtain 

n’<H=nxH^+nxH®  (1.2) 

For  a TM  plane  wave  illumination  the  incident  field  may  be  expressed  as 

5*  - 

5*  - ♦ (1.3) 

and 

nxH^  = is  - HQCos<fe-J^®"°®<^  ^ (l.A) 

From  Eq.  (1.4)  it  is  seen  that  if  we  assume  to  be  independent  of  (j) 
then  H at  the  surface  must  have  a (()  dependence.  However,  H may  be  com- 
puted as  the  field  having  as  its  source.  By  symmetry  a ({>  Independent 
J can  produce  only  a if>  independent  H , and  thus  a contradiction  is  reached. 
It  must,  therefore,  be  concluded  that  the  surface  current  induced  on  a 
nonzero  radius  cylinder  by  an  incident  plane  wave  must  have  a ij>  variation. 
In  addition  to  the  above  analysis,  it  is  Instructive  to  examine  the  solu- 
tion for  an  infinitely  long  cylindrical  scatterer.  It  is  well  known  that 
the  current  distribution  on  the  infinite  long  cylinder  may  be  expressed 
as  a Fourier  series  in  the  angle  variable  if)  [4]. 

The  necessity  for  the  Induced  current  to  have  a (J)  variation  is  often 
overlooked  in  both  theoretical  and  experimental  - ^rk.  The  extent  to 
which  the  angular  variation  affects  the  results  of  course  depends  on  the 
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specific  problem  being  studied.  The  investigation  reported  in  the  follow- 
ing sections  was  initiated  at  this  laboratory  as  a result  of  observations 
made  during  a study  of  a V-cross  scattering  structure  having  ka  = 0.04. 

It  was  observed  that  a 180°  rotation  of  the  structure  and  hence  of  the 
measuring  probe  resulted  in  significant  differences  in  the  measured  cur- 
rents. The  measured  current  distribution  of  the  V-cross  scatterer  are 

shown  in  Figs.  2a  and  2b.  From  this  data  it  can  be  seen  that  the  varia- 

o o 

tion  in  current  between  (J>  = 0 and  ifi  = 180  is  quite  significant  on  some 
portions  of  the  structure.  In  order  to  obtain  a better  understanding  of 
these  effects  a series  of  measurements  on  straight  cylindrical  scatterers 
of  ka  = 0.04  has  been  conducted  with  emphasis  on  determining  the  nature  of 
the  (ji  variation.  The  results  of  this  investigation  are  reported  in  the 
following  sections. 

(ii)  Measurement  Techniques 

The  cylindrical  scatterer  system  used  for  this  investigation  is 
shown  in  Fig.  3.  The  scatterer  is  constructed  from  brass  tubing  of  radius 
3.175mm  mounted  perpendicular  to  a large  vertical  image  plane  (5.5m  by 
8.5m).  The  tube  is  slotted  axially  to  provide  for  movement  of  the  mea- 
suring probe  except  for  the  last  1 cm  at  the  end  of  the  tube.  The  er.d 
of  the  tube  is  unslotted  to  provide  a continuous  conducting  surface  for 
a ^ directed  current  which  may  exist  near  the  end  of  the  tube.  This  test 
structure  is  illuminated  by  a transmitting  antenna  composed  of  a dipole 
mounted  in  a 90°  corner  reflector  located  at  a distance  of  about  10  wave- 
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lengths  from  the  test  structure.  A frequency  of  600 MHz  was  chosen  for 
the  incident  field  giving  a free  space  wavelength  of  X = 50  cm.  At  this 
frequency  the  image  plane  is  IIX  by  17X  and  the  electrical  radius  of  the 
scatterer  is  ka  = 0.04.  The  test  structure  was  mounted  on  a rotatable 
disk  in  the  image  plane.  Rotation  of  the  entire  structure  thus  effec- 
tively places  the  probe  at  different  ^ positions  around  the  cylinder. 

The  current  distribution  was  measured  using  a center  loaded  loop 
probe  3.175mm  diameter  and  constructed  from  Cable-Wave  UT-20  coaxial 

cable.  The  probe  is  oriented  to  measure  H.  and  hence  the  axial  component 

<(> 

of  current.  The  signal  cable  from  the  probe  passes  through  a 1.59mm 
radius  tube  which  is  fastened  to  the  probe  carriage  and  runs  inside  of 
the  test  cylinder  to  instrumentation  behind  the  image  plane.  This  small 
tube  is  coated  with  an  absorbing  material  to  reduce  the  possibility  of 
fields  being  excited  on  the  interior  of  the  test  cylinder. 

During  the  measurements,  the  incident  field  is  monitored  by  a short 
monopole  located  on  the  image  plane  approximately  Im  from  the  test  struc- 
ture and  at  the  same  distance  from  the  transmitter  as  the  test  structure. 
The  signal  received  from  this  probe  is  used  to  correct  the  measurements 
on  both  amplitude  and  phase  for  drifts  in  the  incident  field  illumination. 

Signal  from  the  current  probe  is  measured  using  a HP  8405A  vector 
volt  meter.  The  magnitude  and  phase  outputs  of  the  vector  volt  meter  are 
measured  with  a HP  3490  digital  volt  meter  which  is  interfaced  to  a HP  9820 
calculator.  All  measurements  are  conducted  under  program  control  by  the 
calculator  and  the  data,  after  corrections  for  drifts  of  incident  field, 
are  stored  in  the  calculator  for  later  analysis  and  plotting. 
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As  in  any  measurement  it  is  important  to  understand  the  response  | 

of  the  measurement  probe  and  its  relation  to  the  quantity  being  measured.  i 

i 

A sufficiently  small  center  loaded  loop  antenna  responds  essentially  to 
the  total  magnetic  flux  over  the  area  of  the  loop  [5].  From  the  boundary 
condition  Eq.  (1.1),  it  is  seen  that  the  total  tangential  magnetic  field 
at  the  surface  of  the  conductor  is  equal  to  the  surface  current  density. 

Since  the  signal  induced  in  the  probe  is  taken  as  being  proportional  to 

the  surface  current  density,  errors  may  be  introduced  uecause  of  the  finite 

size  of  the  probe.  The  incident  field  component  is  a very  slowly  varying  ' 

function  over  the  dimensions  of  the  probe.  It  is,  therefore,  only  the 

scattered  components  of  the  total  field  which  may  produce  significant 

errors.  Rapid  spatial  variations  of  the  scattered  field,  however,  can 

only  occur  near  discontinuities  or  changes  in  geometry  of  the  conducting 

surfaces  which  in  the  present  problem  only  occur  near  the  end  of  the  tube. 

Since  measurements  were  made  onlv  to  within  4.25  tube  radii  and  8.5  probe 
radii  of  the  end,  the  errors  over  the  measurement  region  are  small. 

On  all  sliding  surfaces,  i.e.,  probe  carriage  and  rotating  disk, 
a conducting  grease  compound  was  used  to  insure  good  electrical  contact. 

Tests  were  conducted  with  these  areas  covered  by  conducting  tape.  No 
significant  change  in  the  measured  signal  was  observed.  The  conducting 
grease  was,  therefore,  assumed  to  provide  adequate  electrical  continuity. 

Measurements  of  the  surface  current  density  were  made  at  1 cm  intervals 
in  z and  45°  intervals  in  ((>  over  the  surface  of  the  cylinder.  The  result- 


ing data  are  presented  in  the  following  section. 
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(ill)  Measured  Data  and  Analysis 

The  magnitude  and  phase  of  the  measured  current  density  distributions 
for  finite  length  conducting  cylinders  of  ka  = 0.04  are  shown  in  Figs.  4 - 
7 as  a function  of  z,  with  (J>  as  a parameter.  Four  different  lengths  of 
cylinder  were  studied  having  h = X/4,  X/2,  3X/4,  and  X.  Although  data 
were  recorded  at  5°  intervals  in  ij>  for  the  full  360°  range,  only  data  for 
♦ = 0°,  90°,  and  180°  are  presented  since  no  unusual  behavior  was  observed 
at  the  other  angles.  Although  it  can  be  seen  that  the  angular  variation 
of  the  current  is  pronounced  for  all  lengths  studied,  the  effects  may  be 
seen  to  differ  in  some  respects  as  a function  of  h.  Further  discussion 
of  this  h dependence  is  contained  in  the  next  section. 

From  the  solution  for  scattering  by  an  infinite  length  cylinder  [4] 
and  studies  of  an  electrically  thick  cylinder  [6],  it  appears  that  a 
Fourier  series  expansion  in  the  variable  <i>  is  an  appropriate  means  of 

) 

analyzing  the  data.  The  surface  current  density  may,  therefore,  be  ex-  i 

pressed  as 

cos())  + cos  n(t>  + ...  (1.5) 

where  only  cos(n(J>)  terms  are  included  since  these  are  the  only  terms  to  be  1 

expected  if  the  geometry  of  Fig.  3 is  maintained  in  the  experimental  system. 

With  the  aid  of  a computer  the  coefficients  of  Eq.  (1.5)  were  computed 
from  the  measured  data.  The  actual  analysis  Included  computation  of  sin(n(|)) 
coefficients  as  a test  of  the  accuracy  of  the  experimental  system.  In  all 
cases,  the  computed  sin(n(j))  coefficients  were  very  small  and  are  believed 
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to  result  as  much  from  noise  as  from  any  real  features  of  the  measured 

current  density.  The  computed  J (z)  coefficients  are  shown  in  Figs.  8, 

n aw 

9,  10,  and  11. 

(iv)  Conclusions 

At  the  present  time  no  theoretical  solution  for  scattering  from 
finite  length  cylinders  of  moderately  thin  size  is  available  for  comparison 
with  the  measured  data  presented  here.  Some  conclusions  and  understanding 
of  the  behavior  of  the  current  distribution  can  be  drawn,  however,  from  the 
measured  data  and  may  serve  as  a guide  to  the  formulation  of  an  approximate 
theoretical  treatment. 

The  general  behavior  of  the  current  density  is  found  to  be  greatly 
dependent  on  the  length  h of  the  cylinder  as  was  pointed  out  above.  This 
height  effect  is  seen  from  a comparison  of  Figs.  4 and  6 with  Figs.  5 and  7. 
Figs.  4 and  6 for  h = X/4  and  3X/4,  respectively,  show  a 4)  variation  in 
both  magnitude  and  phase  of  the  current  density  but  show  very  similar  z 
dependence  for  all  values  of  <J).  In  contrast,  Figs.  5 and  7 for  h = A/2 
and  A,  respectively,  show  a significant  change  in  z dependence  for  differ- 
ent (|i  positions  especially  near  the  end  of  the  cylinder.  The  most  striking 
feature  of  these  current  distributions  is  that  the  axially  directed  cur- 
rent does  not  approach  zero  toward  the  end  of  the  cylinder.  In  fact,  the 
data  for  both  h = A/2  and  A shows  that  the  magnitude  of  at  iji  = 180° 
reaches  a minimum  and  then  increases  toward  the  end.  The  explanation  for 
the  end  behavior  of  and  the  differences  in  behavior  as  a function  of  h 


is  most  easily  understood  from  an  examination  of  the  Fourier  components 
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shovm  In  Figs.  8 - 11.  For  all  cases  studied  the  J^(z) , i.e.  , cos<ji, 
coefficient  is  essentially  constant  as  is  J2(z).  The  cylinders  of  length 
h = X/4  and  3A/A  are  of  resonant  length  and  the  Jq(z),  i.e.,  angularly 
independent,  component  is  seen  to  be  the  typical  resonant  current  dis- 
tribution which  is  approximated  by  sin  k(h-z) , thus  for  the  resonant  length 
the  JqCz)  component  is  not  only  large  in  magnitude  but  approaches  zero 
toward  the  end  with  a large  slope.  Over  the  range  of  z for  which  data 
are  available  remains  large  compared  to  J^.  The  lengths  h = X/2  and 
X,  on  the  other  hand,  are  antiresonant  and  JqCz)  is  seen  to  be  essentially 
the  forced  response  of  the  general  form  (1  - cos  k(h-z)).  This  forced 
current  distribution  is  not  only  smaller  in  magnitude  relative  to  than 
the  resonant  currents  but  approaches  zero  toward  z = h with  zero  slope. 
Thus,  for  the  antiresonant  cases  becomes  the  largest  component  over  a 
distance  of  about  O.IX  near  the  end.  This  predominance  of  the  component 
near  the  end  explains  the  behavior  of  the  total  current  density  observed 
in  Figs.  5 and  7.  Since  it  is  believed  that  the  and  all  higher  order 
coefficients  will  remain  essentailly  constant  to  z = h,  a sudden  change 
in  the  phase  of  the  current  would  also  be  expected  ver^'  near  the  end  for 
the  resonant  cases,  h = X/4  and  3X/4,  if  it  were  possible  to  obtain  mea- 
surements in  this  region. 

The  constancy  with  respect  to  z of  Jj^(z)  and  higher  order  coefficients 
is  to  be  expected  since  these  components  may  be  considered  as  directly 


excited  by  the  Incident  field  and  must  exist  in  order  to  satisfy  the 
boundary  conditions.  The  behavior  of  these  components  at  the  end  of  the 
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cylinder,  however,  deserves  some  consideration.  From  the  analysis  of 
electrically  thick  cyliider  [6],  it  is  found  that  a (j>  directed  component 
of  current  exists  near  z = h.  It  is  reasonable  to  assume,  therefore,  that 
some  (f)  directed  component  of  current  will  exsit  near  the  end  of  any  non- 
zero radius  cylindrical  scatterer.  Furthermore,  in  accordance  with  the 
edge  condition,  this  (ji  directed  component  must  be  singular  at  z = h. 

Thus,  it  appears  that  the  angularly  dependent  Fourier  components  of  the 
z directed  current  are  essentially  constant  for  all  z;  the  total  current, 
however,  becoming  <()  directed  at  z = h. 

Further  Information  concerning  the  observed  distribution  may  be 
obtained  from  a comparison  with  the  known  solution  for  infinitely  long 
cylinders.  The  induced  current  density  on  a infinite  cylinder  illuminated 
by  a plane  wave  may  be  expressed  as  [4] 


J 

z 


4)  = 


upra 


.-n  ini 
.1  e-" 


(1.6) 


(2) 

where  (ka)  is  the  Hankel  function  of  the  second  kind.  Combining  terms, 
Eq.  (1.6)  becomes; 


2E 

J (,(,)  = — 2- 
z'  mpira 


which  is  of  the  same  form  as  the  Fourier  series  expansion  of  Eq.  (1.5). 
Using  small  argument  approximations  for  the  Ilankel  functions,  the  relative 
magnitude  of  the  Fourier  coefficients  for  ka  = 0.04  are  found  to  be 
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Jq  = 3.38 

j“  = 0.02  j“  (1.8) 

Since  no  absolute  value  of  the  incident  field  illumination  was  available 
for  the  measured  data,  it  is  not  possible  to  compare  the  actual  magnitudes 
of  current  with  the  equivalent  results  for  infinite  cylinder.  However, 
taking  the  values  measured  for  it  is  possible  to  compute  an  equivalent 

OO  00 

value  of  J^,  The  computed  values  of  Jq,  obtained  from  the  average  measured 
vr.iues  of  J^,  are  shown  as  horizontal  dashed  lines  in  Fig.  8 - 11.  It  is 

OO 

seen  that  for  the  antiresonant  lengths  h = X/2  and  X the  equivalent 

is  very  nearly  the  mean  value  of  the  measured  Jq(2)  component.  As  expected, 

00 

however,  the  equivalent  value  of  is  much  less  than  the  mean  value  of 
the  measured  for  the  resonant  lengths  h = X/4  and  3X/4. 

The  results  of  this  investigation  clearly  establish  the  existence  of 
an  angular  variation  of  the  induced  current  on  electrically  thin  cylindrical 
scatterers.  It  is  evident,  therefore,  that,  for  at  least  some  applications, 
solutions  based  on  the  conventional  thin  wire  approximation  will  be  inade- 
quate. Furthermore,  it  is  apparent  that  care  must  be  taken  in  performing 
measurements  on  thin  cylindrical  scatterers.  Measurements  of  the  current 
distribution  which  are  made  at  a single  arbitrarily  chosen  angular  position 
may  be  misleading,  particularly  when  compared  with  thin  wire  solutions. 

The  measured  data,  however,  indicate  that  the  and  higher  order  components 
are  very  small.  Measurements  conducted  at  <{>  = 90°  will,  therefore,  provide 
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sufficient  information  about  the  angularly  independent  component,  J^, 
and  hence  the  total  current.  For  more  complete  information,  the  current 
density  should  be  measured  at  ()>  = 0°  and  ([)  = 180°. 

In  conclusion,  the  observations  drawn  from  the  measured  data  suggest 


that  an  approximate  solution  for  thin  cylindrical  scatterers  might  be 
obtained  by  a combination  of  a thin  wire  approximation  and  the  solution 
for  an  Infinite  length  cylinder. 
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Figure  2b  Phase  of  current  on  V-Cross  Scatterer 


Figure  4a  Magnitude  of  Cu 


Figure  4b  Phase  of  Current  Density 


Figure  5a  Magnitude  of  Current  Density, 
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Figure  6b^  Phase  of  Current  Density,  h = 3X/4 


Figure  7a  Magnitude  of  Current  Density, 


Figure  7b  Phase  of  Current  Density, 


Figure  8 Fourier  Components  of  Current  Density, 
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Figure  10  Fourier  Components  of  Current  Density,  h = 3A/4 
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Figure  11  Fourier  Compom 
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